Flux-Line Lattice Structures in Untwinned YBa 2 Cu 3 07_j 
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A small angle neutron scattering study of the flux-line lattice in a large single crystal of un- 
twinned YBa2Cu307_5 is presented. In fields parallel to the c-axis, diffraction spots are observed 
corresponding to four orientations of a hexagonal lattice, distorted by the a-b anisotropy. A value 
for the anisotropy, the penetration depth ratio, of A a /A(,=1.18(2) was obtained. The high quality of 
the data is such that second order diffraction is observed, indicating a well ordered FLL. With the 
field at 33° to c a field dependent re-orientation of the lattice is observed around 3T. 

PACS numbers: 74.60.Ge, 74.72.Bk, 61.12.Ex 



The remarkable properties of the mixed state in the 
cuprate high-T c superconductors are of great current in- 
terest. In particular, the expectation that high-T c su- 
perconductors have an unconventional pairing symmetry 
has led inevitably to the question: how does the structure 
of the flux-line lattice (FLL) differ between conventional 
and unconventional superconductors? The question has 
been taken up by several recent theoretical contributions 
which predict a variety of interesting FLL effects 
all deviating from the benchmark triangular Abrikosov 
lattice. However, such discussions may presuppose that 
the crystallographic properties of the FLL are already 
well understood in the simplest low-held regime where 
unconventional effects are least prevalent. This has been 
far from the experimental truth. Observations require a 
probe sensitive to the microscopic arrangement of flux- 
lines. Direct imaging |s| and decoration techniques || 
all have inherent drawbacks, and muon spin rotation has 
not yet achieved the sophistication to resolve the most 
subtle effects M. By comparison, small angle neutron 
scattering (SANS) provides unrivalled insights into the 
crystallography of the FLL, and is the only technique 
capable of unequivocably resolving such questions. 

The demanding nature of neutron experiments requires 
large single crystals (masses > 200mg), and because of 
this YBa2Cu307_a (YBCO) has been the cuprate of 
choice for SANS experiments |^|. However, the mate- 
rials properties of YBCO are complicated. The pres- 
ence of Cu-0 chains which are aligned with the crys- 
tallographic b direction render the otherwise tetragonal 
structure orthorhombic. Upon cooling from the growth, 
twin boundaries form along {110} directions separating 
domains of interchanged a and b axes. A strong interac- 
tion between flux-lines and these twin planes significantly 
influences FLL properties. There is a further effect of 
the chains. Although the orthorhombic distortion is only 
slight (w 1%), the electronic structure is markedly af- 
fected, and the consequence is anisotropy within the ab- 



plane of both superconducting (^JT^] and normal state 
PI properties. All previous SANS studies have been on 
heavily twinned crystals ]L^ , |l3| ], and although observa- 
tions of a pattern with four-fold symmetry were claimed 
to be due to unconventional d x 2_ y 2 pairing @, it could 
not be discounted that alignment by twin planes, in com- 
bination with the a-b anisotropy, did not provide a more 
plausible explanation 

A new development in materials technology, involv- 
ing the application of uniaxial stress during the cooling 
process, has made possible the growth of very large un- 
twinned crystals JT^ j . In this paper we report SANS mea- 
surements on such untwinned single crystals of YBCO. 
The results show the effects of the a-b anisotropy on the 
FLL, and prove explicitly that the results from twinned 
crystals [[l2| cannot be interpreted as evidence for d-wave 
effects. Having clarified the low-field picture, we then 
present the first indication at higher fields that the range 
of validity of simple London scaling theory is limited. 

The structure of the lattice formed by flux-lines in a 
superconductor is dictated in crude terms, as a function 
of flux-line density, by the intrinsic anisotropy of the elec- 
tronic structure and the shape of the core of individual 
flux-lines, as well as by the extrinsic effects of flux-line 
pinning by defects and surfaces. The penetration depth, 
A, and coherence length, £, being the parameters which 
define the two length scales of flux-line interactions, their 
ratio k = A/£ is naturally an important quantity. Where 
K is not much bigger than unity, core interactions exist 
over the whole (H,T) space and structures ranging from 
triangular to square are observed, as in borocar bides |l6| ] 
and Nb |L7| , as a manifestation of anisotropic electronic 
structures dependent upon the orientation of the applied 
field with respect to the atomic lattice. In high-K su- 
perconductors, such as cuprates and low-T c NbSe2 Jl8| , 
one would expect that core effects only become influen- 
tial at fields comparable with H C 2, and at low enough 
fields a purely electrodynamic interaction model should 
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be sufficient. In a basic London model, incorporating 
anisotropy in the effective mass tensor, predictions for 
fields applied parallel to the principal c-axis are a hexag- 
onal FLL which is degenerate in energy with respect to 
orientation with the atomic lattice pTgj ] . In reality, the 
presence of even vanishingly small higher-order effects, 
intrinsic or extrinsic in origin, will be able to lift this de- 
generacy to produce a preferred orientation. Among in- 
trinsic effects, it is in the structure of the vortex core that 
an unconventional pairing symmetry is revealed, calcula- 
tions predicting it to be four- fold |§. Free energy calcu- 
lations have therefore been made incorporating four-fold 
symmetry through higher gradient terms in a GL-type 
theory and in a complete microscopic derivation of 
non-local electrodynamics in a London model The 
predictions are for a distorted hexagonal lattice at low 
field transforming to a square lattice over a field range 
which is dependent upon adjustable parameters related 
to the strength of the d x 2_ y 2 contribution. 

The measurements have been obtained over a series 
of experiments using three SANS instruments (Dll, 
D17 and D22) at the Institut Laue-Langevin, Grenoble, 
France. A field was applied parallel to the incident neu- 
tron beam, and the FLL formed by cooling through T c . 
Rotation of the sample about a vertical axis with re- 
spect to the field, and of the field and sample as one 
with respect to the beam were both possible. As in 
all static imaging techniques, we see the FLL structures 
"frozen in" at the irreversibility line near T c . Two crys- 
tals were initially studied, of masses 312mg and 1125mg, 
grown and detwinned under the same conditions, and 
oxygenated in flowing O2 at 490°C for optimal doping. 
Magnetization measurements showed T c of 92K. Essen- 
tially identical results were obtained from both, and the 
work presented therefore relates solely to the larger. 

In Figure |l], a diffraction pattern is shown, obtained 
at 1.5K with a relatively low field of 0.51T applied par- 
allel to the c-axis. The diffraction spots lie in an ellip- 
tical ring, representing a FLL which is essentially poly- 
crystalline. This is in distinct contrast to patterns ob- 
tained for twinned crystals [ ^2|Jl^1 which show four-fold 
patterns. The ellipse shape is clearly due to the super- 
conducting anisotropy between the a and b basal plane 
directions. For this to be observable requires that the 
bulk of the sample must have a single Cu-0 chain orien- 
tation. For fields larger than B c i ^0|, the axial ratio of 
the ellipse is equal to the ratio of the penetration depths, 
Jab = A Q /A& (making = m a /rrib the effective mass 
ratio). The sign of the anisotropy is consistent with a 
reduced penetration depth Af, due to the increased su- 
percurrent flow along the chain direction. From fits to 
the ellipse of scattering, a value for 7 of 1.18(2) was ob- 
tained. This is slightly higher than the 1.13 estimated 
by Bitter decoration |(| , but in the lower range of values 
found by polarized reflectivity, 7=1.3-1.6 (9), and Joseph- 
son tunnel junctions, 7=1.2-1.7 [jlO|. That the variability 
is between samples rather than techniques is borne out 
by /zSR measurements on the same sample which give 
7=1.16(2) (2l]]. No field dependence in 7 was observed. 



Reflectivity measurements on strongly oxygen disordered 
samples give 7 as low as 1.05 p^| . We conclude that the 
Cu-0 chains in our crystal are slightly disrupted by a 
small concentration of oxygen vacancies or impurities. 

The intensity distribution around the ring in Fig. ^, 
although unbroken, contains distinct diffraction spots. 
From the assessment of many such patterns, we conclude 
that four different flux-line lattice orientations are being 
observed, each contributing a hexagonal pattern of six 
spots distorted by the a- b anisotropy. The situation is 
illustrated in Fig Two of the structures have a FLL 
plane oriented with the axes of the atomic lattice, and 
two have a plane at 45° to them (that is with the {110} 
directions). Surprisingly, the latter two are the same as 
two of the four orientations observed in twinned crys- 
tals. The implication is that this crystal is not com- 
pletely free of twins, although the domains of the op- 
posite orientation must be a very small fraction of the 
volume otherwise a second set of spots aligned on an el- 
lipse at right angles would be observed. A neutron study 
of the crystallography of these YBCO crystals has con- 
firmed this ]23|] , the minority orientation was found to 
be confined to less than 5% of the sample. The two 
45° distorted hexagonal patterns correspond to part of 
the square pattern observed by Keimer et al |l2] |, who 
claimed that the orientation and distortion were due to 
d-wave effects. It is made clear by the present observa- 
tions that the FLL structures, instead, are controlled at 
these low fields by anisotropy plus twin plane pinning in 
the manner postulated earlier JlJ] . Although a non-local 
contribution to the penetration depth anisotropy is pos- 
sible Ipij , the field- dependence of the flux latice distortion 
argues against this being important. 

By tilting the c-axis ~ 1.5° (about the vertical) away 
from the applied field direction the different nature of the 
FLLs is made apparent. As expected of flux- lines pinned 
along their length by extended defects, the two 45° FLLs 
remain fixed to the twin plane direction while the re- 
maining two lie instead along the field direction. Con- 
sequently in Fig. ||, this remarkably clear pattern shows 
only the twelve first order spots of the lattices oriented 
with the a and b axes. A ring of second order reflections 
is testament that the FLL is well formed. The left-hand 
side is stronger because these spots have been brought 
closer to their Bragg condition by rotating the field and 
sample as one. The Bragg condition for spots of the 45° 
FLLs are separated a further 1.5° away, which is consid- 
erably larger than the rocking curve widths. Gaussians 
fits to individual spots gave full-widths at half-maximum 
of 0.65(5)° and 0.71(5)° for the a-b and 45° FLLs respec- 
tively (instrumental broadening contributes only 0.2°). 
These give a lower limit for the longitudinal correlation 
lengths of the flux- lines of 1.8/im. 

The intensity of a reflection integrated over its full 
rocking curve can be related to the penetration depth. 
The established expression for the structure factor de- 
rived from London theory |25| indicates that the intensity 
is proportional to A~ 4 . Making this calculation for the 
rocking curve of the pattern in Fig. □ requires that we 
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estimate the proportion of the total number of flux-lines 
contained in each lattice. The rocking curves give (within 
10%) equal proportions of the intensity distributed be- 
tween them. Using this, we obtain A=138(5)nm. If we 
take A= \/\ a Xb, and use the already measured value of 
7, A a =150(6)nm and Ab=127(6)nm are obtained. The 
value of A a provides a useful check of the consistency of 
our results, as it is known to be insensitive to the state 
of the chain layers, and is therefore reproducible across 
samples. It is typically quoted as 155nm for optimal dop- 
ing |26f| which agrees within experimental error. 

Over the range of fields investigated, up to a maximum 
of 4 Tesla, the FLL struture remained essentially un- 
changed. This may not be inconsistent with the various 
unconventional and non-local theories whose distortions, 
if present, may simply be too weak over this range. We 
note that the distorted square lattices observed by STM 
at higher fields S are not consistent with our observa- 
tions. Our first indication of the breakdown of simple 
London scaling theory comes from measurements made 
with the field inclined at an angle to the c-axis. 

Rotating the field away from c to large angles, another 
factor becomes important: the anisotropy between su- 
percurrents along c compared with those in the a6-plane. 
Previous SANS measurements on twinned crystals have 
given 7 C = \ c /\ a b ~4.5 |b|. This j c anisotropy mixes 
with the jab anisotropy. With b as the vertical axis of 
rotation the eccentricity of the ellipse is accentuated by 
7 C adding to j a b, and at angles comparable to 30° from 
c a continuous highly eccentric ellipse of scattering is 
observed, indicating a polycrystalline FLL. In contrast, 
when a is the axis of rotation, spots from the two ab 
aligned FLLs remain, with the a-oriented FLL becoming 
dominant with increasing angle. This is the orientation 
expected from simple London theory [jl9| . A possible ex- 
planation for the difference rotating about the two axes, 
is that rotation about b leaves the chains perpendicular 
to the applied field, while rotation about a leads to a shal- 
lower angle between field and chain direction, suggesting 
an unexpected influence of the chain structure upon the 
flux-lines. At an angle of «33° from c, shown in Fig. 
^(a), the two anisotropies cancel, and an almost undis- 
torted hexagonal pattern is observed, again with second 
order spots. The FLL is a-axis oriented; however, under 
the same conditions but at a field of 3.8T, the pattern 
in Fig. |](b) is observed corresponding to a FLL oriented 
instead with the 6-axis. There is clearly a field depen- 
dence to whatever mechanism is controlling the FLL ori- 
entation, or perhaps a crossover from the London mech- 
anism at low flux-line density to another at high density. 
At present we have no explanation for the effect, and it 
may yet prove necessary to invoke core effects, perhaps 
of d x 2_ y 2 origin, to account for it. 

In conclusion, we emphasise that the crystallography 
of the FLL must be considered in the context of a compli- 
cated variety of influences, before addressing the singular 
question of unconventional pairing. The results demon- 
strate the various implications for the FLL caused by 
the presence of CuO chains in YBCO; introducing mi- 



crostructural complications of twinning, a basal plane 
anisotropy in the effective mass, and even the suggestion 
of an interaction between the flux-lines and the chain 
structure. To resolve whether the observed field depen- 
dant reorienation is an intrinsic feature of unconventional 
superconductivity will require experiments on cuprates 
free of chains and ideally with tetragonal symmetry. 
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FIG. 1. The SANS pattern with 0.51T applied parallel to 
the c-axis; a-axis is vertical in the figure. The plot area cov- 
ers 128x128 detector pixels, in reciprocal space T0.0173A -1 ; 
neutron wavelength X n — 12A, detector distance 14.5m, and 
collimation 14.4m. 
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FIG. 2. Illustration of how the pattern in Fig. Ill is formed 
from two reciprocal lattices having planes oriented with the 
primary atomic lattices (top), and two having planes aligned 
with the {110} and {110} twin planes (bottom) producing 
spots at 45° (emphasised). 



FIG. 3. The pattern with a field of 0.2T, not quite parallel 
to c. A ring of twelve first-order spots from two hexagonal 
lattices is visible, and remarkably well defined second-order 
spots. The plot covers 64x64 detector pixels, ±0.0144A _1 in 
reciprocal space; A n = 10A, detector distance 14m, and colli- 
mation 13.6m. 
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FIG. 4. The applied field tilted towards the fe-axis at an 
angle of 33° to the c-axis. At 0.2T in the upper pattern a 
single hexagonal lattice aligned with the a-axis (vertical in 
the plots) is observed; ±0.0144A , A„ = 10A, detector dis- 
tance 14m, and collimation 13.6m. At fields above 3T, the 
lower pattern changes to be oriented instead with the 6-axis; 
±0.0503A , A n =10A, detector distance and collimation are 
4 2.5m. 



